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1.  Introduction 


1.1.  Background 

Pneumatic  artificial  muscle  actuators  were  first  developed  in  artificial  limb 
research  in  the  1950’s  and  1960’s  (Gavrilovic  and  Marie,  1969,  Schulte,  1961),  They 
were  commercialized  by  the  Bridgestone  Corporation,  Japan,  in  the  1980’s  (Inoue,  1988) 
using  rubber  as  the  primary  muscle  material  for  robotics  applications  (e.g.,  Feng,  et  al, 
1997,  Hesselroth,  et  al,  1994,  Uno,  1986,  Yamafuji,  et  al,  1993).  These  actuators  have 
high  power-to-weight  and  power-to-volume  ratios  and  flexible  mechanical  features 
friendly  to  human  as  compared  to  hydraulic  and  electric  actuators.  Therefore,  these 
actuators  have  shown  considerable  potential  for  enhancing  human  strength  and  mobility 
for  military  as  well  as  civilian  applications. 

However,  some  features  of  the  pneumatic  actuators,  namely  the  compressibility  of 
the  air  inside  the  muscles,  the  delay  of  the  response  and  its  nonlinearity  for  example, 
make  it  difficult  to  precisely  control  output  displacements  and  forces  of  these  actuators, 
which  prevents  their  effective  applications  in  the  fields  as  mentioned  above. 

Significant  efforts  have  recently  been  made  to  overcome  such  drawbacks.  For 
example,  Repperger,  et  al,  (1997)  performed  static  and  dynamic  tests  on  a  pneumatic 
muscle  actuator  system  to  identify  its  nonlinear  characteristics  and  to  incorporate  them 
into  a  force  control  algorithm.  Chou  et  al  (1996)  developed  a  linearized  model  of  the 
McKibben  artificial  muscle  pneumatic  actuator.  Galdwell  et  al  (1995)  applied  the 
adaptive  control  scheme  combined  on  on-line  system  identification  to  improve  the 
accuracy  of  the  position  control.  Hesselroth,  et  al  (1994)  developed  a  neural  learning 
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algorithm  for  position  control  of  a  robot  arm  actuated  by  rubber  muscles.  The  previous, 
work  represented  by  these  papers  contributed  significantly  to  the  artificial  muscle 
technology.  Additional  work,  however,  is  needed  to  further  improve  the  control 
performance  of  pneumatic  muscles  in  order  to  make  them  implementable. 


1.2.  Objective  of  this  Research 

The  objectives  of  this  research  are 

i)  to  develop  a  prototype  pneumatic  muscle  actuator  system  reflecting  the  state  of 
the  art  technology. 

ii)  to  study  and  model  its  static  and  dynamic  characteristics  and 

iii)  to  develop  and  apply  effective  control  algorithms  to  achieve  accurate  position 
control. 
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2.  Development  of  a  Pneumatic  Muscle  Actuator  System 
2. 1 .  Characteristics  of  Pneumatic  Air  Muscles 

The  pneumatic  muscles  used  as  actuators  in  this  research  are  called  McKibben 
Pneumatic  Artificial  Air  Muscles.  The  muscle  has  an  internal  rubber  bladder  covered 
with  an  external  shell  formed  from  interwoven  nylon  fibers  (Figure  2.1).  When  the 
internal  bladder  is  pressurized,  the  high  pressure  air  pushes  against  its  inner  surface  and 
its  external  shell,  and  thus  tends  to  increase  its  volume.  The  bladder  and  thread  are  fixed 
at  both  ends  tightly  with  fittings.  Due  to  the  non  extensibility  (or  very  high  longitudinal 
stiffness)  of  the  nylon  fibers  of  the  external  shell,  the  muscle  shortens  according  to  its 
volume  increase  and/or  tension  force  if  it  is  coupled  to  a  mechanical  load.  By  combining 
a  pair  of  such  muscles  in  different  ways  and  by  controlling  the  pressure  difference  of  the 
air  inside  the  two  muscle  bladders,  one  can  provide  various  actuations. 

Praldad  Threads 


Figure  2.1:  Structure  of  McKibben  Pneumatic  Artificial  Muscle 

Three  types  of  pneumatic  muscles  were  used  for  the  experiment.  The  specifications 
of  each  muscle  are  shown  in  Table.  2.1. 
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Table  2.1:  Specifications  of  each  Muscle 


2.2.  Experimental  Setup  of  One-Linked  Arm  Fixture 


An  experimental  setup  of  a  one-linked  arm  fixture  was  designed  and  fabricated. 
As  shown  in  Figure.  2.2,  the  system  consists  of  a  pair  of  muscles,  an  air  compressor,  two 
proportional  pressure  control  valves,  an  air  regulator,  an  encoder,  two  pressure  sensors,  a 
rotary  link,  and  a  computer  with  AD/DA  boards.  The  muscles  are  connected  through  a 
pulley,  which  represents  an  arm  joint.  The  angular  position  of  the  joint  is  controlled  by 
adjusting  the  longitudinal  deformation  of  each  air  muscle  through  controlling  the  air 
pressure  inside  the  muscle.  The  encoder  is  used  to  measure  the  angular  position  of  the 
joint  and  the  two  pressure  sensors  are  used  to  monitor  the  air  pressure  inside  each  muscle. 
The  measured  data  from  these  sensors  are  sent  to  the  computer  through  the  analog  to 
digital  (A/D)  converter.  The  pressure  control  signals  are  sent  to  the  air  valve  from  the 
computer  via  the  digital  to  analog  (D/A)  converter.  Figure  2.3  is  a  photo  of  the  fixture. 


Figure  2.2;  Schematic  Diagram  of  Fixture 
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Figure  2.3:  Photo  of  Fixture 


2.3.  Specifications  of  Devices  used  in  the  Experiment 

The  specifications  of  the  valves,  pressure  sensors,  encoder,  A/D  and  D/A 
converters  and  the  air  compressor  are  shown  in  Table  2.2,  through  Table  2.6.  Figures  2.4 
through  2.8  show  the  photos  of  these  devices. 
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Table  2.2:  Specifications  of  Magnetic  Proportional  Pressure  Control  Valve 


|  Model 

PPC5B-AAA-AGCA 

Manufacturer 

Fluidtrols  Corporation 

Pressure  Range 

100  (psi) 

Pressure  Reference 

Gage  Pressure 

Overall  Accuracy 

□  0.5%  F.S. 

0-10  V  Differential 

Figure  2.4:  Pressure  Valve 
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Table  2.3:  Specifications  of  Pressure  Sensor 


Model 

DP2-42N  | 

Manufacturer 

SUNX 

Type  of  Pressure 

Gage  Pressure 

Rotate  Pressure  Range 

0  ~  l.OOOMpa 

Applicable  Fluid 

Non-corrosive  gas 

Supply  Voltage 

Output  Voltage 

1  to  5  V 

Zero-point:  within 

IV  □  5%  F.S. 

Span:  within  4V  0  5%  F.S. 
Linearity:  within  V  []  5%  F.S. 

Figure  2.5:  Pressure  Sensor 


Table  2.4:  Specifications  of  Rotary  Encoder 


Model 

L2SG-F 1-2500-P4C- 
8830-LED-SM16 

Manufacturer 

BEI  Sensors  & 

Systems  Company 

Code 

Incremental 

10,000 

Maximum  RPM 

EEEIiEHHHHHH 

Moment  of  Inertia 

EmUESSESHi 

Voltage  Supply 

+5V  de 

2  channels  (cw  &  ccw) 

bbhi 

100  kHz 

i’*?" '  •  ‘T’ ?  *  r-  '■  ■■ 


Figure  2.6:  Rotary  Encoder 
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Table  2.5:  Specifications  of  Multifunction  I/O  Board 


1  Model 

PCI-6024E 

Manufacturer 

National  Instruments 

Analog  Input 

16  Single-ended  or  8 
differential 

200kS/s  sampling  rate 
200kS/s  stream-to-disk 
rate 

12  bit  resolution 

Analog  Output 

2channels 

12-bit  resolution 

8  (5V/TTL)  lines 

Counter  /  Timer 

2  up/down 

24  bit  resolution 

Triggering 

Digital 

Figure  2.7:  Multifunction  I/O  Board 


Table  2.6  Specifications  of  Air  Compressor 


Model 

WL340003AJ 

Manufacturer 

CAMPBELL 

HAUSFELD 

Maximum  Pressure 

Container  Volume 

Figure  2.8:  Air  Compressor 
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3.  Characteristics  of  Experimental  Fixture 
3. 1.  Theoretical  Modeling  of  Air  Muscle 

In  order  to  control  the  angular  position  or  the  torque  of  the  arm  joint  through 
controlling  the  pressure  of  the  air  muscle,  the  relationship  between  the  muscle  force,  air 
pressure  and  the  contraction  ratio  of  the  muscle  must  be  identified.  In  this  research,  a 
theoretical  relationship  was  first  derived.  In  order  to  simplify  the  model  the  following 
assumptions  were  made. 

•  The  shape  of  the  braided  sleeve  is  cylindrical. 

•  The  threads  in  the  braided  mesh  shell  are  not  extensible. 


•  The  elasticity  of  bladder  is  ignored. 

•  The  friction  of  each  thread  is  ignored. 

As  depicted  in  Figure  3.1,  the  length  of  the  muscle,  L,  has  the  following  relationships 
with  the  contraction  ratio  e,  and  the  initial  length  L„. 


L  =  (\-s)Ln 
L'  =  (l-s)L'n 


(3.1) 


where 

L  =  Length  of  the  muscle 
L„  =  Initial  length  of  the  muscle 

L'  =  Longitudinal  length  of  the  muscle  for  one  turn  of  thread. 

L'  -  Initial  longitudinal  length  of  the  muscle  for  one  turn  of  thread. 
As  such 


Ln  —  mLn 
L  =  ml! 

where  m  =  number  of  turns  of  thread 


(3.2) 
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Figure  3.1:  Simple  Model  of  Air  Muscle 


The  relationship  between  the  above  parameters  are  illustrated  in  the  triangles 


shown  in  the  Figure  3.2 


The  geometry  of  the  parameters  of  the  triangle  is  as  follows 

(nDn)2+L'„2  =(nD)2+L'2  =  l2  (3.3) 

where  Dn  and  D  are  at  initial  and  final  conditions  of  the  air  muscle,  respectively. 

By  solving  for  D,  yields 

D2  =  A.2  +~V{l-(l-eJ2}  (3-4) 

7t 

Using  the  same  geometry  and  solving  for  the  circumference,  yields 

KDn=L'ntm<j)n  (3  5) 

nD  =  L'tan<j> 

where,  (j)n  and  <|>  refer  to  the  gradient  angle  of  each  thread  at  initial  and  contracted  states. 
Thus,  setting  the  circumference  equal  to  each  other  and  solving,  yields 


D  (1  -e)Dn 
taiuj)  tan</>„ 

(3.6) 

Similarly, 

Ln  =  ml  cos  <j>n 

(3.7) 

and 

c 

k  -35 

II 

(3.8) 

which  yields. 

L  nD„ 

tan  f, 

(3.9) 

By  substituting  Equation  3.9  into  Equation  3.4,  yields, 

1  tan2^  J 

(3.10) 
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where, 

Afc  =  the  circular  force  that  act  on  a  small  portion  of  the  air  muscle. 

P  =  the  pressure  inside  the  air  muscle. 

r,  D  =  the  radius  and  the  diameter  of  the  muscle  respectively. 


Assuming  Fc  is  the  circular  Force  that  is  generated  in  the  length  of  a  one  turn  of  thread, 
then 


Fc  =  A fcL' 


Thus, 


2  Fc  =  PDL' 


(3.16) 


Thus, 


jr.  PDL' 

Fr  = - 

c  2 

The  longitudinal  force  generated  by  internal  air  pressure  Flp  is  shown  as: 


(3.17) 


Flp=-EP  (3.18) 

where,  the  force  pointing  towards  the  inside  of  the  muscle  is  assumed  the  positive 
direction,  and 

E  =  cross  section  of  the  air  muscle. 


E  = 


(3.19) 


Thus, 

Flr=~nD2P  (3.20) 

and,  the  contraction  force  of  air  muscle,  F  is  the  balancing  force  as  shown  in  Figure  3.4. 
F  =  Ft+F„  (3.21) 
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(3.22) 


F,  =  F  +  -nD2P 


Now  substituting  Equation  3.17  and  Equation  3.22  into  Equation  3.14,  yields 


r^p.rnL 

4  2tan0 


(3.23) 


Substituting  Equation  3.12  into  Equation  3.23  yields  the  theoretical  Equation, 


F  +  -kD2P  =  —  -f-^- 
4  2  l^tan</> 


(3.24) 


By  substituting  Equation  3.24  into  Equation  3.6,  yields 


*P,2Q-g)2 

2  tan 2  (j>n 


F  +  —nD2 P  -  '  ~  J  P 


(3.25) 


By  substituting  Equation  3.25  into  Equation  3.10,  yields 


f  =  ^TT-{3cos20,(l-e)!-l} 
4sm  <j>n 


(3.26) 


Thus,  the  theoretical  model  of  the  relationship  of  the  muscle  force,  contraction  ratio  and 
the  pressure  can  be  expressed  as 


F  =  A(B(l-s)2-l)P 


(3.27) 


where, 


4  sin 2  (f>n 
B  =  3  cos2  <f>n 


(3.28) 
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3.2.  Experimental  Modeling  of  Air  Muscle 

The  theoretical  model  in  Equation  3.27,  which  shows  the  relationship  between  the 
force,  pressure  and  contraction  ratio  of  each  muscle  is  obtained  in  the  Section  3.1. 
However,  theoretical  and  experimental  data  do  not  agree  perfectly  with  each  other  as 
shown  in  Figures.  3.5  a,  b,  and  c.  The  two  curves  of  the  experimental  results  are  for  the 
forward  and  backward  motions  respectively.  The  force  value  for  each  contraction  ratio 
from  the  experiments  is  generally  lower  than  the  value  from  the  theoretical  model.  The 
difference  may  be  due  to  such  reasons  as  the  assumption  of  uniform  cross  section  of  the 
muscle,  the  friction  being  ignored  for  each  thread,  and  the  elasticity  of  the  bladder  being 
ignored. 

Therefore  the  experimental  models  are  also  obtained  for  each  of  the  muscles, 
RUB-125,  AM-02  and  AM-04.  Figures  3.6,  3.8,  and  3.10  show  the  best-fitted  straight 
lines  of  the  experimental  data  for  three  different  pressure  values. 

As  shown  in  Figure  3.6,  the  experimental  data  approximates  a  straight  line  within 
the  tested  range  of  the  pressure  and  contraction  ratio  (shown  in  Table  3.1)  when 
hysteresis  is  ignored.  Furthermore,  the  slope  of  the  line  is  proportional  to  the  pressure. 
The  force  value  at  e  =  0  is  also  in  proportion  to  the  pressure. 


Table  3.1:  Condition  for  Experiment 


1 

RUB- 125 

AM-02 

AM-04 

Pressure  Range  [psi] 

20-40 

20-40 

20-50 

Range  of  Contraction  Ratio  s 

0.02-0.196 

0.167-0.233 

0.20  -  0.270 

Total  Stroke  Range  of  8 

0.0 -0.2 

0.0  -  0.3 

0.0  -  0.3 
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Therefore,  the  experimental  model  showing  the  relationship  among  the  muscle 
force,  contraction  ratio  and  pressure  can  be  expressed  in  the  Equation  shown  below. 

F  =  (aP  +  b)e+(cP  +  d )  (3.29) 

where  P  =  pressure 

a,b,c,  and  d  are  constants. 

Table  3.1  lists  the  contraction  and  pressure  ranges  for  each  type  of  muscles  for 
which  the  force-contraction  ratio  relationship  can  be  linearized.  Figures  3.7,  3.9  and  3.11 
show  the  amount  of  error  between  the  above  mentioned  experimental  data  and  their  best 
fit  linear  lines. 

For  the  RUB- 125  muscle,  Figure  3.6  shows  the  linear  relationships  at  three 
pressure  states;  28.4,  42.6,  56.8[psi]  (  2,  3,  4  [kgf/cm2]  ).  The  maximum  stroke 
(contraction  ratio)  of  the  RUB-125  muscle  is  0.2.  The  range  of  the  contraction  ratio 
tested  in  the  experiment  is  0.02  to  0.196,  which  covers  97  [%]  of  the  entire  range. 
Choosing  the  parameters,  a,  b,  c  and  d  in  the  Equation  3.29  to  best  fit  the  experimental 
data,  the  error  between  the  approximated  line  and  the  experimental  data  of  force  values 
can  be  minimized  below  15  [%],  as  shown  in  Figure.  3.7. 
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Figure  3.5  (a)  &  (b):  Relationship  between  Force  and  Contraction  Ratio 
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:  Error  between  Model  and  Experimental  Data  (RUB-125) 
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Figure  3.8:  Experimental  Modeling  (AM-02) 
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Figure  3.9:  Error  between  Model  and  Experimental  Data  (AM-02) 
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Figure  3.11:  Error  between  Model  and  Experimental  Data  (AM-04) 
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3.3.  Parameter  Identification  of  Air  Muscle 


The  identification  of  the  parameters  of  air  muscles,  were  done  experimentally 
using  the  fixture  shown  in  chapter  2.  The  simple  depiction  of  the  system  is  shown  in 
Figure.  3.12. 


Figure  3.12:  Construction  of  Air  Muscle  Actuator  System 

In  this  discussion  the  displacements  were  assumed  to  be  small.  As  mentioned 
previously,  the  Equation  of  force  vs.  contraction  for  each  i  th  muscle  is  given  by 

F,  =  ( aP,  +  b)s,  +(cPt+d)  (3.30) 

3.3.1.  The  Equation  of  Motion 

With  the  assumption  that  Pj  and  8j  change  infinite  decimal  amounts  from  the  basic 
pressure  Pq  and  basic  contraction  ratio  s0  (e  under  P  =  Po), 


Pt=P»+8P, 

(3.31) 

=  £0f  +  i 

(3.32) 
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(3.33) 


o  *  r6 
0£j  =  —OS  2  = - 

Z/„ 


where, 

R  =  The  radius  of  the  pulley 
0  =  The  angle  of  rotation 
Ln  =  The  natural  length  of  the  air  muscle 
Substituting  these  Equations  into  Equation  3.30  yields, 

F,  =  (aPofiot  +  beo,  +  cP0l  +  d)  +  (ae0i  +  c)8P,  +  (aP0i  +  b)8ei  +  a8Pi8ei 
If  it  is  assumed  that 
F^Fu+SF, 

then, 

8Ft  =  Ka-8Pi  +  Kb-8sl 

where, 

F0i  —  ^Poi^oi  +  bs0j  +  cP0l  +  d 

Ka  =  ae0i  +  c 
Kb  =  aP0j  +  b 


(3.34) 


(3.35) 


(3.36) 


(3.37) 


abPfiSt  «0 

Using  these  Equations,  the  dynamic  Equation  of  motion  can  be  obtained  as  follows. 

The  moment  M  at  the  pulley  with  the  link  is 

M  =  (Fl-F2)-r  (3.38) 

where,  Fj  and  F2  are  the  contraction  force  of  muscle  1  and  muscle  2  respectively. 

If  the  force  of  each  muscle  under  basic  pressure  is  F01  =  F02  ,then 
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M  =  (8Fl  —8F2)r 


(3.39) 


Substituting  Equation  3.36  into  Equation  3.39  yields 
M  =  {Ka(8Pl  -8P2)  +  Kb(8ei  - 8s 2 )}  •  r 
Substituting  Equation  3.33  into  Equation  3.40  yields 


M  =  Ka-r(8P]  -8P2)  +  Kb-r 


r2  rO^ 


V  Ln  J 


thus 


7  Kh 

M  =  Ka-r(8Pl -5P2)+  6 

Finally,  the  moment  M  of  the  rotating  parts  can  be  described  as 
M  =  Kc(8Px  -SP2)~  Kd  -9 

where, 

Kc  =  Ka-r 

=  - 

L, 

The  Equation  of  motion  is  shown  as 
Jd+cB=M 

where, 

J=  The  moment  of  inertia  of  the  wheel  and  the  link 
C=  Damping  coefficient 


Substituting  Equation  3.43  into  Equation  3.45  yields 
J9  +  c6  +  KdQ  =  Kc(8Px  -  8P2 ) 


(3.40) 


(3.41) 


(3.42) 


(3.43) 


(3.44) 


(3.45) 


(3.46) 
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3.3.2.  Parameter  Identification  by  Experiment 

The  parameters  Kd  and  c  in  Equation  3.46  can  be  obtained  by  using  the  free 
vibration  of  our  system.  The  system  was  subjected  to  an  initial  angle  while  keeping  the 
pressure  of  both  muscles  at  basic  pressure  Po.  The  back  and  forward  motion  of  the  link 
was  recorded.  The  results  of  each  muscle  are  shown  in  Figure. 3.1 3. 

If  both  muscles  are  kept  at  the  same  constant  pressure,  the  Equation  3.46  can  be 
expressed  as 


JG+cd  +  KdO  =0 
Dividing  the  Equation  by  J  yields 

6  +2<Zcon9  +(o„26  =0 


where, 


C;.  Damping  ratio 
con:  Natural  frequency 


Kd  =  Jco„2 
c  =  2  J<X 


Also 


©«  = 


2  K 


Substituting  these  into  Equation  3.49  and  Equation  3.50  respectively,  yields 


Kd  = 


Ak2J 

T2 


c  = 


4 <  J 
T 2 


(3.47) 

(3.48) 


(3.49) 

(3.50) 

(3.51) 

(3.52) 

(3.53) 
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Tim  e  [s  ] 
RUB-125  muscle 


0,4  0,6  0. 

Tim  e  [s  ] 

AM-02  muscle 


2  0.4  0.6  0 

Tim  e  [s] 


AM-04  muscle 

Figure  3*13:  Free  Vibration  of  Link 


The  period  T  and  the  damping  ratio  C,.  are  obtained  using  the  recorded  data  of 
free  vibration.  Then  substituting  these  parameters  into  Equation  3.52  and  Equation  3.53 
the  parameters  Kd  and  c  are  obtained. 

The  parameter  Kc  was  determined  by  performing  the  static  tests  to  obtain  the 
relationship  between  the  difference  of  pressure  change  (5Pi  -8P2)  of  each  muscle  and  the 
angle  of  0. 


If  the  moment  M  =  0  then  the  Equation  3.43  can  be  simplified  as 

Kc=-™-e 

8PX  -  SP2 

if  P01  =  P02  in  Equation  3.31  then 
A P  =  8Pl-8P2=Pl-P2 
Therefore,  Equation  3.54  becomes 

Kc  =  —  •  Kd 
A P 


(3.54) 


(3.55) 


(3.56) 


The  experimental  results  obtained  from  the  relationship  between  0  and  AP,  are 
shown  in  Figure.  3.14.  The  values  of  parameters  obtained  from  three  types  of  muscles 
using  experimental  results  are  shown  in  Table  3.2. 


Table  3.2:  Parameter  Values  Identified  by  Free  Vibration  Experiment 


Kd[Nm/rad]  | 

c  [Nm/(rad/s)] 

RUB- 125 

82.2 

EEKEi^H 

AM-02 

12.7 

KBHI9H 

AM-04 

114.5 

0.605 

-1.77  10‘6 
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4.  System  Control  and  Response 


With  the  analytical  model  established  and  its  parameters  experimentally  identified, 
algorithms  were  developed  for  controlling  the  angular  displacement  of  the  link  though 
position  control  of  a  pair  of  pneumatic  muscles.  First,  Proportional  Integral  (PI)  control 
was  applied  which  resulted  in  a  time  delay  of  desired  response.  Then  a  Feed-Forward 
control  algorithm  was  developed  and  applied,  which  successfully  compensated  the  time 
delay  in  muscle  response. 

4.1.  PI  Control 

The  encoder  measured  the  angular  displacement  of  the  pulley  and  this  signal  was 
sent  to  the  computer.  The  computer  then  sent  a  pressure  control  signal  to  the  proportional 
pressure  control  valves  of  the  two  muscles  based  on  the  difference  of  the  measured  angle 
and  a  reference  angle.  The  control  signal  is  based  on  the  PI  control  shown  as  follows. 

8P„  =  KPQe  +  K,  feedt  (4.1) 

Prl  =/>01  +8P  /2 

(4.2) 

Pn=P01-5Pr/2 

where, 

6e  —6r  —Q:  difference  between  the  reference  angle  and  the  measured  angle 

Pri,  Pr2*  pressure  control  signals  for  muscle  1  and  muscle  2 
Kp :  proportional  gain 
Ki :  integral  gain 

The  output  voltage  signals  sent  to  the  pressure  control  valves  are 
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r,  =  Kvprl 


where, 

Ky  (i=l,  2):  coefficient  of  the  valve  based  on  the  relationship  between  input 
voltage  and  output  pressure 

Step  responses  under  the  PI  control  were  obtained  through  experiments  and 
compared  with  the  results  simulated  using  the  experimental  model  shown  in  Equation 
3.29.  Figure  4.1  show  the  experimental  and  simulated  results  for  each  type  of  the 
muscles,  where  the  reference  signal  was  set  at  5  degrees.  The  parameter  values  used  in 
the  experiments  and  simulation  are  shown  in  Table  4.1. 

In  the  case  of  RUB- 125  and  AM-02  muscles,  the  experimental  and  simulated 
responses  agreed  well,  except  for  their  frequencies.  Ignoring  the  friction  of  threads  and 
the  elasticity  of  bladder  when  modeling  the  muscle  actuator  system,  mainly  caused  the 
difference  in  frequencies. 

In  case  of  the  largest  size  muscle,  AM-04,  the  experimental  result  did  not  well 
correlate  with  the  simulated  one.  Due  to  the  large  air  capacity  of  the  muscle  and  the 
small  inlet  to  supply  air,  the  response  velocity  of  muscle  was  limited,  which  was,  not 
considered  in  the  analytical  model.  From  the  experimental  data,  it  can  be  confirmed  that 
the  maximum  contraction  velocity  under  the  basic  pressure  around  35[psi]  is  8.64[mm/s]. 
Therefore,  for  a  step  control  signal,  the  required  response  velocity  exceeded  the 
maximum  velocity  capacity  of  the  muscle  and,  as  a  result,  the  muscle  contracts  at  its 
maximum  velocity. 

In  summary,  the  simulated  response  reasonably  agreed  well  with  the  experimental 
response,  indicating  that  the  analytical  model  together  with  its  parameters  identified  in 
this  study  is  sufficient. 
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Table  4.1:  Parameter  Values  for  PI  Control 


PARAMETERS  ||  RUB-125  ““ |  AM-02  |  AM-04  ] 

J  [kgm']  “ 

0.029 

r[m] 

0.045 

Ln  [m] 

0.160 

0.358 

C  [Nm/(rad/s)] 

0.476 

0.243 

0.605 

T  [s] 

0.118 

0.3 

0.1 

C 

0.154 

0.2 

0.166 

Kc[m3] 

IjgQQI 

Kd  [Nm/rad] 

82.22 

12.72 

114.49 

Kp[[psi/rad] 

65 

120 

100 

Kj[[psi/rad.sec] 

1000 

1300 

800 
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4.2.  Compensation  using  Feed-Forward  Control 


This  section  documents  the  results  of  feed-forward  control  that  was  used  to 
control  the  movement  of  the  link  arm  in  a  predetermined  path.  The  path  sought  was  a  sine 
curve.  Time  delay  of  air  muscle  system  was  shown  in  Figure  4.1.  A  feed-forward  control 
algorithm  was  thus  developed  to  compensate  the  time  delay. 

The  control  system  equipped  with  the  feed-forward  algorithm  is  shown  in  Figure. 
4.2  (Noritsugu  et  al,  1986).  The  gain  Kfp  was  decided  by  the  static  relationship  between 
the  angular  displacement  and  the  difference  of  pressures  AP,  which  was  discussed  in 
Chapter  3.3.2.  In  addition,  a  low  pass  filter  was  added  with  a  time  constant  T  =  0.125 
seconds  to  alleviate  the  rapid  change  in  velocity  of  motion.  Adding  a  differential  motion 
with  gain  of  Kfv  also  helped  to  compensate  the  delay. 
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The  angular  displacement  of  the  link  arm  was  controlled  to  follow  a  sinusoidal 
curve  with  an  amplitude  of  5  degrees  and  a  period  of  2  seconds  for  muscles  RUB- 125  and 
AM-02,  and  4  seconds  for  muscle  AM-04.  The  Figure  4.3  shows  the  response  of  each 
type  of  muscle  system  under  the  PI  control  without  feed-forward  compensation.  In  this 
case  when  compared  with  the  reference  signal,  time  delay  in  response  is  again 
demonstrated. 

The  sinusoidal  response  under  the  feed-forward  control  is  shown  in  Figure,4.4. 
The  parameter  values  of  the  control  system  are  shown  in  Table  4.2.  Considering  the 
limited  velocity  capacity  of  muscle  AM-04,  sinusoidal  signals  with  higher  frequencies 
(period  less  than  4  seconds)  were  not  tested.  Compared  with  the  results  under  the  PI 
control  without  feed-forward  compensation  shown  in  Figure.  4.3,  the  results  under  the 
feed-forward  control  were  significantly  better  in  terms  of  the  reduction  of  time  delay. 
The  maximum  error  between  the  reference  and  measured  signals  is  reduced  to  about  0.3 
degrees.  The  instability  demonstrated  at  the  early  stages  of  the  time  history  is  due  to  the 
discontinuity  of  the  velocity  at  the  starting  point. 

Table  4.2:  Parameter  Values  for  Control  with  Feed-Forward  Compensation 


1 

RUB-125 

AM-02 

AM-04 

Kfp  [psi/radl 

100 

245 

102 

Kfy  fpsi*s/radl 

20 

45 

45 

Kp  [psi/radl 

100 

120 

1500 

Kj  [psi/(rad*s)l 

20 

280 

3800 
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Figure  4.3:  Experimental  Results  of  PI  Control 


5.  Summary  and  Conclusions 


The  following  tasks  have  been  performed  in  this  research  project: 

1.  A  prototype  pneumatic  rubber  muscle  actuator  system  has  been  designed  and 
fabricated  which  consists  a  pair  of  muscles  and  a  pulley,  resembling  a  human  arm 
joint.  An  encoder  was  installed  to  measure  the  rotation  of  the  pulley  driven  by  the 
muscles  and  two  pressure  sensors  to  measure  the  pressure  inside  the  rubber 
blades.  Air  was  supplied  from  an  air  compressor  through  two  servo  values.  This 
system  was  connected  to  a  computer  through  I/O  boards  including  A/D  and  D/A 
converters. 

2.  Three  types  of  rubber  muscles  with  different  manufacturers  and  different  sizes 
have  been  tested. 

3.  The  theoretical  model  of  the  muscle  system  (relationship  among  the  pressure, 
contraction  ratio  and  force)  was  found  quite  different  from  the  experimental  data. 
As  a  result,  experiments  were  performed  to  develop  a  model  and  to  identify  its 
parameter  values,  rather  than  using  the  theoretical  model.  The  developed  model 
was  used  for  designing  the  control  system. 

4.  Control  algorithms  were  developed  including  a  feed-forward  algorithm  to 
compensate  the  time  delay  of  the  muscle  system. 

The  following  observations  and  conclusions  can  be  made  from  this  research: 

1 .  The  linearity  of  the  muscle  actuator  depends  on  the  manufacturing  quality.  Based 
on  the  criterion  determined  in  this  study,  the  linear  range  of  the  Bridgestone 
muscles  was  found  to  be  as  large  as  97  [%]  of  its  stroke,  while  the  that  of  the  two 
Shadow  muscles  only  reached  22  [%]  and  23  [%]  of  their  strokes  respectively. 
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2.  The  accuracy  of  the  analytical  model  of  the  muscle  system  experimentally 
developed  in  this  study  was  demonstrated  through  comparison  of  the  experimental 
and  analytical  step  responses. 

3.  The  time  delay  of  the  muscle  response  depends  on  the  size  of  the  muscle.  The 
larger  muscle  suffers  from  longer  time  delay,  although  it  provides  larger  power 
and  stroke. 

4.  The  control  algorithm  equipped  with  a  feed-forward  compensation  developed  in 
this  research  is  very  effective  in  reducing  the  time  delay  in  response. 

5.  Excellent  control  performance  of  the  pneumatic  muscle  actuator  system  was 
demonstrated  by  experiments. 
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